Wideband transceiver for antenna array

ABSTRACT

Beam measurement values are computed by transforming a signal to a beamspace representation using a beamforming matrix that includes a plurality of columns that define a plurality of beams. A second plurality of beams is selected from the plurality of beams based on a power captured by each beam derived from the computed beam measurement values. For each beam of the second plurality of beams as a first beam, a spatial phase shift component and a temporal delay component is computed based on the fixed spatial angle associated with the first beam and a carrier frequency, a filter is defined using the spatial phase shift component and the temporal delay component, a beam measurement value associated with the first beam from the computed beam measurement values is selected, and the defined filter is applied to the selected beam measurement value to define a filtered value.

REFERENCE TO GOVERNMENT RIGHTS

This invention was made with government support under 1247583 awarded by the National Science Foundation. The government has certain rights in the invention.

BACKGROUND

Wideband high-dimensional antenna arrays are expected to play a key role in future wireless systems. Due to narrow beamwidths, phased array/beamforming methods are the natural choice for design and analysis of high-dimensional multiple-input, multiple-output (MIMO) systems. However, these methods are based on a narrowband assumption that is violated as the bandwidth and array dimension increase.

SUMMARY

A method of processing a communication signal is provided. A plurality of beam measurement values are computed by transforming a signal received from each antenna of a plurality of antennas to a beamspace representation using a beamforming matrix that includes a plurality of columns that define a plurality of beams. Each column of the plurality of columns defines a steering vector for a beam that corresponds to a fixed spatial angle relative to a boresight vector extending from the plurality of antennas. A second plurality of beams is selected from the plurality of beams based on a power captured by each beam of the plurality of beams derived from the computed plurality of beam measurement values. For each beam of the selected second plurality of beams as a first beam, a spatial phase shift component is computed based on the fixed spatial angle associated with the first beam and a carrier frequency, a temporal delay component is computed based on the fixed spatial angle associated with the first beam and the carrier frequency, a filter is defined using the spatial phase shift component and the temporal delay component, a beam measurement value associated with the first beam is selected from the computed plurality of beam measurement values, and the defined filter is applied to the selected beam measurement value to define a filtered value. The defined filtered values are combined.

In another example embodiment, a receiver is provided that includes a processor configured to perform the method of processing a communication signal.

In yet another example embodiment, a transmitter is provided that includes a plurality of antennas configured to receive a first signal and to transmit a second signal and a processor configured to compute a plurality of beam measurement values by transforming the first signal to a beamspace representation using a beamforming matrix. The beamforming matrix includes a plurality of columns, wherein each column of the plurality of columns defines a steering vector for a beam that corresponds to a fixed spatial angle relative to a boresight vector extending from the plurality of antennas. The plurality of columns thereby define a plurality of beams. The processor is further configured to select a second plurality of beams from the defined plurality of beams based on a power captured by each beam of the plurality of beams that is derived from the computed plurality of beam measurement values. For each beam of the selected second plurality of beams as a first beam, the processor is further configured to compute a spatial phase shift component based on the fixed spatial angle associated with the first beam and a carrier frequency, to compute a temporal delay component based on the fixed spatial angle associated with the first beam and the carrier frequency, and to define a filter using the computed spatial phase shift component and the computed temporal delay component. The processor is further configured to filter a transmit signal through each defined filter, to process the filtered signals through the beamforming matrix, and to transmit the processed, filtered signals from the plurality of antennas as the second signal.

Other principal features of the disclosed subject matter will become apparent to those skilled in the art upon review of the following drawings, the detailed description, and the appended claims.

BRIEF DESCRIPTION OF THE DRAWINGS

Illustrative embodiments of the disclosed subject matter will hereafter be described referring to the accompanying drawings, wherein like numerals denote like elements.

FIG. 1 depicts a communication system in accordance with an illustrative embodiment.

FIG. 2 depicts a transmitting transceiver and a receiver array of a receiving transceiver of the communication system of FIG. 1 in accordance with an illustrative embodiment.

FIG. 3 depicts a block diagram of the receiving transceiver of FIG. 1 in accordance with an illustrative embodiment.

FIG. 4 depicts a beamspace of the receiving transceiver of FIG. 1 in accordance with an illustrative embodiment.

FIGS. 5a-5e show plots of S_(i)(f) for an M=61 antenna array for i_(o)=25 as a function of a normalized frequency for five values of i centered on i_(o). FIG. 5a shows S_(i)(f) for i=25; FIG. 5b shows S_(i)(f) for i=24; FIG. 5c shows S_(i)(f) for i=26; FIG. 5d shows S_(i)(f) for i=23; and FIG. 5d shows S_(i)(f) for i=27.

FIG. 6a shows a plot of Σ_(i) S_(i)(f) for the M=61 element array for three beams centered at i_(o), i={24, 25, 26}.

FIG. 6b shows a plot of Σ_(i) S_(i)(f) for the M=61 element array for five beams centered at i_(o), i={23, 24, 25, 26, 27}.

FIG. 7 shows an optimal receiver implementation of the receiving transceiver in accordance with an illustrative embodiment.

FIG. 8 shows a near optimal, low-complexity receiver implementation of the receiving transceiver in accordance with an illustrative embodiment.

FIG. 9a presents signal to interference and noise ratio (SINR) curves as a function of signal to noise ratio (SNR) for a phased array receiver with different values for a channel dispersion factor in accordance with an illustrative embodiment.

FIG. 9b presents SINR curves as a function of SNR for a three-beam receiver with the different values for the channel dispersion factor in accordance with an illustrative embodiment.

FIG. 9c presents SINR curves as a function of SNR for a five-beam receiver with the different values for the channel dispersion factor in accordance with an illustrative embodiment.

FIG. 10 shows a near optimal, low-complexity transmitter implementation of the transmitting transceiver in accordance with an illustrative embodiment.

DETAILED DESCRIPTION

Referring to FIG. 1, in an illustrative communication system, there is a line-of-sight (LoS) path between a first transceiver 100 and a second transceiver 102 that represents clear spatial channel characteristics though first transceiver 100 and second transceiver 102 also may be linked in a multipath environment. For example, a signal 104 transmitted by second transceiver 102 is radiated towards first transceiver 100 on the LoS path. First transceiver 100 and second transceiver 102 support both the transmission and the reception of electromagnetic waves. Use of the terms transmitter and receiver is to describe an example function that can be performed by each device. For purposes of discussion, second transceiver 102 is denoted as a transmitting transceiver, and first transceiver 100 is denoted as a receiving transceiver though each transceiver may be configured to support either or both functions. First transceiver 100 is illustrated as a base station of a communications system and second transceiver 102 is illustrated as a communications device that communicates with the base station such as a cell phone.

One or both of first transceiver 100 and second transceiver 102 may be mounted on moving objects such that a distance between the transceivers may change with time. As known to a person of skill in the art, the communication environment between first transceiver 100 and second transceiver 102 may fluctuate due to changes in environmental conditions such as weather, due to changes in interference sources, and due to movement between first transceiver 100 and second transceiver 102, which may change the multipath environment, any of which may cause a fluctuation in the received signal-to-noise ratio (SNR), signal-to-interference ratio (SIR), and/or signal to interference and noise ratio (SINR) even where the transmission power and other signal characteristics such as frequency, pulsewidth, bandwidth, etc. remain unchanged.

Second transceiver 102 may include a single antenna or a plurality of antennas arranged to form an array. Referring to FIG. 2, first transceiver 100 may include a plurality of antennas 200 arranged to form an array 202. For illustration, array 202 may be a uniform or a non-uniform linear array, a rectangular array, a circular array, a conformal array, etc. The plurality of antennas 200 are mounted in a common plane. An antenna of the plurality of antennas 200 may be a dipole antenna, a monopole antenna, a helical antenna, a microstrip antenna, a patch antenna, a fractal antenna, a feed horn, a slot antenna, etc. An antenna spacing 204, which may be denoted d, may separate each of the plurality of antennas 200 from an adjacent antenna of the plurality of antennas 200 in the common plane. The plurality of antennas 200 are configured to receive an analog signal from second transceiver 102 and/or to radiate a plurality of radio waves toward second transceiver 102. The plurality of antennas 200 may include a number of antennas that may denoted as M.

A boresight vector 206 extends from a center of array 202 perpendicular to the common plane in which the plurality of antennas 200 is mounted. Second transceiver 102 is located along a direction vector 208 which defines an angle 210, which may be denoted φ_(o), relative to boresight vector 206. For illustration, φ_(o) represents only the azimuth angle relative to a linear array. Alternative embodiments can be extended to two-dimensional arrays in which the angle φ_(o) is replaced by a pair of angles representing the azimuth angle and the elevation angle.

Referring to FIG. 3, a block diagram of a receiver 300 of first transceiver 100 is shown in accordance with an illustrative embodiment. Second transceiver 102 may include similar elements as understood by a person of skill in the art. Receiver 300 may include the plurality of antennas 200, a processor 302, a computer-readable medium 304, and a beamforming matrix 306. Different and additional components may be incorporated into receiver 300.

Processor 302 may form a signal that includes one or more information symbols that is sent to second transceiver 102. Processor 302 further may receive a received signal that includes one or more information symbols from the plurality of antennas 200. The received signal may be received by the plurality of antennas 200 after transmission by second transceiver 102. Processor 302 may implement a variety of well-known processing methods, collectively called space-time coding techniques, which can be used for encoding information into information symbols. Processor 302 further may perform one or more of converting a data stream from an analog to a digital form and vice versa, encoding the information symbols, decoding the information symbols, modulating the information symbols, demodulating the information symbols, up-converting the information symbols from a carrier frequency, that may be denoted f_(c), down-converting the information symbols to the carrier frequency, performing error detection and/or data compression, Fourier transforming the information symbols, inverse Fourier transforming the information symbols, etc.

Processor 302 may be implemented as a special purpose computer, by logic circuits, by hardware circuits, in firmware, and/or using any combination of these methods. Processor 302 further may execute instructions meaning it performs/controls the operations called for by that instruction as understood by those skilled in the art. The term “execution” is the process of carrying out the operation called for by an instruction. The instructions may be written using one or more programming language, scripting language, assembly language, etc. Processor 302 operably couples with computer-readable medium 304. Processor 302 may retrieve a set of instructions from a permanent memory device and copy the instructions in an executable form to a temporary memory device that is generally some form of random access memory (RAM). Receiver 300 may include a plurality of processors that use the same or a different processing technology.

Computer-readable medium 304 is an electronic holding place or storage for information so the information can be accessed by processor 302 as understood by those skilled in the art. Computer-readable medium 304 can include, but is not limited to, any type of RAM, any type of read only memory (ROM), any type of flash memory, etc. such as magnetic storage devices (e.g., hard disk, floppy disk, magnetic strips, . . . ), optical disks (e.g., compact disc (CD), digital versatile disc (DVD), . . . ), smart cards, flash memory devices, etc. Receiver 300 may include one or more computer-readable media that use the same or a different memory media technology. For example, computer-readable medium 304 may include different types of computer-readable media that may be organized hierarchically to provide efficient access to the data stored therein as understood by a person of skill in the art. As an example, a cache may be implemented in a smaller, faster memory that stores copies of data from the most frequently/recently accessed main memory locations to reduce an access latency.

Referring to FIG. 4, first transceiver 100 may be configured to communicate with a plurality of second transceivers 400 that may be of the same or different types. For example, the plurality of second transceivers 400 may include various devices that transmit signals wirelessly using any type of standardized or proprietary communication protocol, frequency, bandwidth, etc. Example devices include computing devices of any form factor including smart phones, tablets, laptops, desktops, servers, etc., vehicles, sensing elements, etc.

The plurality of second transceivers 400 may be distributed over a region 402 randomly, uniformly, non-uniformly, etc. For example, cell phone devices may be distributed randomly throughout region 402. In the illustrative embodiment of FIG. 4, a non-uniform distribution of the plurality of second transceivers 400 is shown. Array 202 of first transceiver 100 is located sufficiently far from region 402 such that far-field assumptions apply. Region 402 can be spatially divided into M spatial beams 404 and a plurality of time delay rings 406 based on the spatio-temporal resolution supported by the plurality of antennas 200 of first transceiver 100. Region 402 may span azimuth angles of ±180 degrees and elevation angles of ±90 degrees.

For illustration, a single-input, multiple-output (SIMO) communication between second transceiver 102 comprised of a single antenna transmitter and first transceiver 100 comprised of an M-dimensional (M receive antennas) uniform linear array (ULA) receiver operating at a carrier frequency f_(c) can be defined. As used herein, x is a vector; X is a matrix; x(t) is a vector valued function of time with Fourier transform X(f)=

{x(t)}=∫x(t)e^(−j2πft)dt. A two-sided signal bandwidth for the SIMO communication can be defined as W=αf _(c),αε(0,2]  (1) where α is a fractional bandwidth and is typically α<<1. The transmitter of second transceiver 102 may transmit a packet of duration T:

$\begin{matrix} {{s(t)} = {\sum\limits_{l = 0}^{N - 1}{s_{l}{\Psi_{l}(t)}}}} & (2) \end{matrix}$ where the {s_(l)}_(t=0) ^(N-1) are independent information symbols with energy E[|s_(l′)|²]=

and the {ψ_(l)(t)}_(l=0) ^(N-1) form an orthonormal basis for the N≈TW dimensional signal space.

In contrast with separable phased array based models, the M-dimensional complex baseband signal r(t) received at array 202 can be related to s(t) transmitted by second transceiver 102 by

$\begin{matrix} {{\underset{\_}{r}(t)} = {{{\left( {{\underset{\_}{h}}^{*}s} \right)(t)} + {\underset{\_}{w}(t)}} = {{\sum\limits_{\ell = 0}^{N - 1}{{s_{l}\left( {{\underset{\_}{h}}^{*}\Psi_{l}} \right)}(t)}} + {\underset{\_}{w}(t)}}}} & (3) \end{matrix}$ where h(τ) is an M×1 spatial channel impulse response and w(t) is a spatially and temporally white complex additive white Gaussian noise with power spectral density N_(o). Taking the Fourier transform yields

$\begin{matrix} {{\underset{\_}{R}(f)} = {{F\left\{ {\underset{\_}{r}(t)} \right\}} = {{\sum\limits_{l = 0}^{N - 1}{s_{l}{\underset{\_}{H}(f)}{\Psi_{l}(f)}}} + {\underset{\_}{W}(f)}}}} & (4) \end{matrix}$ where H(f)={

{h(τ)} is a spatial frequency response.

The beamspace channel representations in time and frequency are obtained by projecting the signal onto a set of orthonormal array steering vectors (M spatial beams 404) at first transceiver 100. A steering vector a_(M)(θ) is defined as

$\begin{matrix} {{a_{M}(\theta)} = \left\lbrack {\mathbb{e}}^{{- {j2\pi}}\; k\;\theta} \right\rbrack_{k \in {I{(M)}}}} & (5) \end{matrix}$ where I(M)=l−(M−1)/2:l=0, . . . , M−1} is a symmetric set of indices centered around 0. Columns of beamforming matrix 306, U_(M), are steering vectors that correspond to M fixed spatial angles, illustrated by M spatial beams 404, with uniform spacing Δθ=1/M.

$\begin{matrix} {U_{M} = {\frac{1}{\sqrt{M}}\left\lbrack {a_{M}\left( {i\;{\Delta\theta}} \right)} \right\rbrack}_{i \in {I{(M)}}}} & (6) \end{matrix}$ and that represent M orthogonal beams (M spatial beams 404) forming a basis for the M-dimensional spatial signal space. Beamforming matrix 306, U_(M), is a unitary discrete Fourier transform (DFT) matrix U_(M) ^(H)U_(M)=U_(M)U_(M) ^(H)=I_(M). The beamspace system models in time and frequency are

$\begin{matrix} {{{\underset{\_}{r}}_{b}(t)} = {{U_{M}^{H}{\underset{\_}{r}(t)}} = {{\sum\limits_{l = 0}^{N - 1}{{s_{l}\left( {{\underset{\_}{h}}_{b}^{*}\Psi_{l}} \right)}(t)}} + {{\underset{\_}{w}}_{b}(t)}}}} & (7) \end{matrix}$

$\begin{matrix} {{{\underset{\_}{R}}_{b}(f)} = {{U_{M}^{H}{\underset{\_}{R}(f)}} = {{\sum\limits_{l = 0}^{N - 1}{s_{l}{{\underset{\_}{H}}_{b}(f)}{\Psi_{l}(f)}}} + {{\underset{\_}{W}}_{b}(f)}}}} & (8) \end{matrix}$ where the beamspace channel impulse response and the beamspace channel frequency response are given by h _(b)(τ)=U _(M) ^(H) h (τ), H _(b)(f)=U _(M) ^(H) H (f).  (9)

The signal s(t) arrives at each antenna (array element) of array 202 with a slightly different delay. Assuming without loss of generality that the delay is 0 at the array center, the delay at the k^(th) antenna of array 202 is

$\begin{matrix} {{\tau_{k} = {k\;{\delta\tau}}},{{k \in {I(M)}};{{\delta\tau} = {\frac{d}{c}{\sin\left( \varphi_{o} \right)}}}}} & (10) \end{matrix}$ where δτ is the delay between adjacent antennas and c is the speed of light. Down mixing the signal s(t) from f_(c), also manifests a delay as a phase shift between array antennas, θ_(k)=kθ_(o), determined by a normalized spatial angle

$\begin{matrix} {\theta_{o} = {{{\delta\tau}\; f_{c}} = {\frac{d}{\lambda_{c}}{{\sin\left( \varphi_{o} \right)}.}}}} & (11) \end{matrix}$ As understood by a person of skill in the art, λ=c/f_(c).

In narrowband systems, the effects of the delay on s(t) can be ignored resulting in the baseband phased array model: r (t)=a _(M)(θ_(o))s(t)+ w (t).  (12) However, as the array dimension and bandwidth increase, this is no longer true and the channel's spatial impulse response is

$\begin{matrix} {{\underset{\_}{h}(\tau)} = \left\lbrack {h_{k}(\tau)} \right\rbrack_{k \in {I{(M)}}}} & (13) \\ {{h_{k}(\tau)} = {{\mathbb{e}}^{{- {j2\pi}}\; k\;\theta_{o}}W\;\sin\;{c\left( {W\left( {\tau - {k\;{\delta\tau}}} \right)} \right)}}} & (14) \end{matrix}$ where bandlimiting the spatial impulse response leads to the sincs. The k^(th) antenna of the spatial frequency response is H _(k)(f)=

{h _(k)(τ)}=e ^(−jπk(δτf+θ) _(o))  (15)

From equations (5), (11), and (15), the channel's spatial frequency response in the general case is H (f)=a _(M)(θ(f)),fε[−W/2,W/2]  (16) with frequency dependent spatial angle

$\begin{matrix} {{\theta(f)} = {{\theta_{o}\left( {\frac{f}{f_{c}} + 1} \right)}.}} & (17) \end{matrix}$

The beamspace channel impulse response is

$\begin{matrix} {{{{\underset{\_}{h}}_{b}{()}} = {{U_{M}^{H}{\underset{\_}{h}(T)}} = \left\lbrack {h_{b,i}(T)} \right\rbrack_{i \in {I{(M)}}}}}{{{h_{b,i}{()}} = {{\frac{1}{\sqrt{M}}{a_{M}^{H}\left( {i\;\Delta\;\theta} \right)}{\underset{\_}{h}(T)}} = {\frac{1}{\sqrt{M}}{\sum\limits_{k \in {I{(M)}}}{{\mathbb{e}}^{{- {j2\pi}}\;{k{({\theta_{O} - {i\;{\Delta\theta}}})}}}W\;\sin\;{c\left( {W\left( {T - {k\;\sigma\; T}} \right)} \right)}}}}}},}} & (18) \end{matrix}$ and the beamspace channel frequency response is

$\begin{matrix} {{{{\underset{\_}{H}}_{b}(f)} = {{U_{M}^{H}{\underset{\_}{H}(f)}} = \left\lbrack {H_{b,i}(f)} \right\rbrack_{i \in {I{(M)}}}}}{{H_{b,i}(f)} = {{\frac{1}{\sqrt{M}}{a_{M}^{H}\left( {i\;\Delta\;\theta} \right)}{\underset{\_}{H}(f)}} = {\frac{1}{\sqrt{M}}{{D_{M}\left( {{\theta(f)} - {i\;\Delta\;\theta}} \right)}.}}}}} & (19) \end{matrix}$

Here

${D_{M}(\theta)} = \frac{\sin\left( {\pi\; M\;\theta} \right)}{\sin({\pi\theta})}$ is the Dirichlet sinc function with D_(M)(0)=M and D_(M)(mΔθ)=0 for integers m≠0.

While the baseband phased array model from equation (12) shows no dispersion in spatial angle or time, equations (14) and (16) show dispersion in time across an aperture of the plurality of antennas 200 and dispersion in spatial angle across the bandwidth. Thus, a point source in narrowband systems spreads across spatial angle and time in wideband, high dimensional systems. Not accounting for the dispersion can lead to severe power loss and pulse distortion as shown below.

The magnitude of the channel dispersion can be determined by how many resolvable beams and delays the channel occupies. The angular spread of the channel over the bandwidth is Δθ_(ch)=|θ(W/2)−θ(−W/2)|=α|θ_(o)|  (20)

Dividing this by the orthogonal beam spacing Δθ=1/M yields Δθ_(ch) /Δθ=Mα|θ _(o)|,  (21) so the channel spans approximately Mα|θ_(o)| orthogonal beams. Similarly, the delay spread across the aperture is Δτ_(ch) =M|δτ|.  (22) Dividing this by delay resolution Δτ=1/W yields Δτ_(ch) /Δτ=WΔτ _(ch) =Mα|θ _(o)|,  (23) so the channel spans approximately Mα|θ_(o)| resolvable delays. Thus, a channel dispersion factor Δ_(ch) =Mα|θ _(o)|=Δτ_(ch)/Δτ=Δθ_(ch)/Δθ  (24) captures the spreading of the channel in both angle and delay. For transmitters located at φ=0 (in the direction of boresight vector 206), from equation (11), θ_(o)=0 so Δ_(ch)=0 and no dispersion occurs. However, for any other angle, θ_(o)≠0, dispersion occurs and increases as the bandwidth and array dimension increase.

The effects of Δ_(ch) can be illustrated by considering an energy spectral density of the channel impulse response:

$\begin{matrix} {{S_{\underset{\_}{h}}(f)} = {{{{\underset{\_}{H}}^{H}(f)}{\underset{\_}{H}(f)}} = {{{{\underset{\_}{H}}_{b}^{H}(f)}{{\underset{\_}{H}}_{b}(f)}} = {{\sum\limits_{i \in {I{(M)}}}{{H_{b,i}(f)}}^{2}} = {{\sum\limits_{i \in {I{(M)}}}{S_{i}(f)}} = M}}}}} & (25) \end{matrix}$ where S_(i)(f)=|H_(b,i)(f)|² is an energy density captured by the i^(th) beam. The channel power is

$\begin{matrix} {\sigma^{2} = {{\frac{1}{W}{\int_{- \frac{W}{2}}^{\frac{W}{2}}{S_{\underset{\_}{h}}(f)}}} = {{\sum\limits_{i \in {I{(M)}}}\sigma_{i}^{2}} = M}}} & (26) \end{matrix}$ where σ_(i) ² is the power captured by the i^(th) beam:

$\begin{matrix} {\sigma_{i}^{2} = {\frac{1}{W}{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{S_{i}(f)}{{\mathbb{d}f}.}}}}} & (27) \end{matrix}$

Because S _(h) (f) is constant with respect to frequency, if a receiver captures all of the channel power there is no pulse distortion. However, if a receiver does not capture all of the channel power, the energy spectral density captured by the receiver may not be flat and pulse distortion occurs. This can be seen by considering the case when θ_(o)=i_(o)Δθ for some i_(o)εI(M). From equations (17) and (19), in this case,

$\begin{matrix} {{S_{i}(f)} = {\frac{1}{M}{D_{M}^{2}\left( {{\theta_{o}\frac{f}{f_{c}}} + {\left( {i_{o} - i} \right){\Delta\theta}}} \right)}}} & (28) \end{matrix}$ and from equation (24) the channel dispersion parameter is Δ_(ch)=αi_(o). In particular,

${S_{i_{o}}(f)} = {{{{a_{M}^{H}\left( \theta_{o} \right)}{\underset{\_}{H}(f)}}}^{2} = {\frac{1}{M}{D_{M}^{2}\left( {\theta_{o}\frac{f}{f_{c}}} \right)}}}$ is the power captured by a phased array receiver.

FIGS. 5a-5e show plots of S_(i)(f) for an M=61 element array for i_(o)=25 (φ_(o)=55°) as a function of normalized frequency for five values of i centered on i_(o). FIG. 5a shows S_(i)(f) for i=25; FIG. 5b shows S_(i)(f) for i=24; FIG. 5c shows S_(i)(f) for i=26; FIG. 5d shows S_(i)(f) for i=23; and FIG. 5d shows S_(i)(f) for i=27. When Δ_(ch)<<1 (|f/f_(c)|<<0.02 over the bandwidth), S_(i) _(o) (f)≈M and S_(i≠i) _(o) (f)≈0. So the i_(o) ^(th) beam has a flat energy density (over the bandwidth of interest), and a phased array receiver exhibits no power loss or pulse distortion as shown in FIG. 5a . Increasing Δ_(ch) causes variation in S_(i) _(o) (f), in particular for Δ_(ch)=2 (|f/f_(c)|≦0.04 over the bandwidth), as shown in FIG. 5a , nulls begin to appear in S_(i) _(o) (f). Thus, using a phased array receiver for larger values of Δ_(ch) results in severe power loss and pulse distortion. On the other hand, from FIGS. 5b -5 e, it is clear that the power lost in S_(i) _(o) (f) is present in the beams adjacent to the i_(o) ^(th) beam.

As shown in FIG. 6a , for Δ_(ch)≦(|f/f_(c)|≦0.06 over the bandwidth), combining the three beams centered on i_(o) (FIG. 5a-5c ) results in an approximately flat power density equal to M (see also equation (25)). Thus, a receiver that uses these three beams results in nearly no power loss or pulse distortion. If Δ_(ch) is increased to 5 (|f/f_(c)|≦0.1 over the bandwidth), the three beams are no longer able to capture all the power, again resulting in pulse distortion and power loss. However, using all five beams shown in FIGS. 5a-5e results in an approximately flat power density as shown in FIG. 6b . In general, for a channel with channel dispersion factor Δ_(ch), a receiver using approximately Δ_(ch) beams centered at i_(o) is able to capture approximately all the channel power and to introduce almost no pulse distortion.

The sufficient statistics {z_(l)}_(l=0) ^(N-1) for detecting the information symbols {s_(l)} are obtained by taking an inner product of the received signal with waveforms {g _(l)(t)}_(l=0) ^(N-1) representing a mapping of

{ψ_(ℓ)(t)}_(ℓ = 0)^(N − 1) into the spatio-temporal receive signal space. The inner product is < x,y>=∫y ^(H)(t) x (t)dt=∫Y ^(H)(f) X (f)df  (29) with the associated norm ∥x∥²=<x, x>. This inner product may be calculated in either the spatial domain or beamspace, and in time or frequency. Thus, the z_(l) are given by

$\begin{matrix} {z_{\ell} = {\left\langle {\underset{\_}{r},{\underset{\_}{g}}_{\ell}} \right\rangle = {{{\sum\limits_{\ell^{\prime} = 0}^{N - 1}{s_{\ell^{\prime}}\left\langle {\left( {{\underset{\_}{h}}^{*}\psi_{\ell^{\prime}}} \right),{\underset{\_}{g}}_{\ell}} \right\rangle}} + \left\langle {\underset{\_}{w},{\underset{\_}{g}}_{\ell}} \right\rangle} = {{s_{\ell}A_{\ell}} + {\sum\limits_{\ell^{\prime} \neq \ell}{s_{\ell^{\prime}}B_{\ell,\ell^{\prime}}}} + {W_{\ell}.}}}}} & (30) \end{matrix}$ where the signal amplitude and interference are A _(l)=<( h*ψ _(l)),g _(l) >,B _(l,l′)=<( h*ψ _(l′))g _(l))  (31) and W_(l)˜

(0,N_(o)∥g_(l)∥²) represents the noise. For a given choice of basis functions {ψ_(l)(t)}, the SINR for the l^(th) test statistic (z_(l)) is

$\begin{matrix} {{SINR}_{\ell} = {\frac{\frac{E_{s}}{N_{o}}{A_{\ell}}^{2}}{{\frac{E_{s}}{N_{o}}{\sum\limits_{\ell^{\prime} = \ell}{B_{\ell,\ell^{\prime}}}^{2}}} + {k_{\ell}}^{2}}.}} & (32) \end{matrix}$

In the low SNR limit

$\left. \frac{E_{s}}{N_{o}}\rightarrow 0 \right.$

$\begin{matrix} {\left. {SINR}_{\ell}\rightarrow{SNR}_{\ell} \right. = {\frac{E_{s}}{N_{o}}\frac{{A_{\ell}}^{2}}{{k_{\ell}}^{2}}}} & (33) \end{matrix}$ which represents the signal power captured by the receiver. Conversely in the high SNR limit

$\left. \frac{E_{s}}{N_{o}}\rightarrow\infty \right.$

$\begin{matrix} {\left. {SINR}_{\ell}\rightarrow{SIR}_{\ell} \right. = \frac{{A_{\ell}}^{2}}{\sum\limits_{\ell^{\prime} \neq \ell}{B_{\ell,\ell^{\prime}}}^{2}}} & (34) \end{matrix}$

A natural choice for g _(l)(t) is a matched filter

$\begin{matrix} {{{\underset{\_}{g}}_{\ell}(t)} = {\frac{1}{\sqrt{M}}\left( {{\underset{\_}{h}}^{*}\psi_{\ell}} \right)(t)}} & (35) \end{matrix}$

with norm ∥g _(l)∥²=1, which results in signal amplitude

$\begin{matrix} {A_{\ell} = {{\frac{1}{\sqrt{M}}{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{\psi_{\ell}^{*}(f)}{a_{M}^{H}\left( {\theta(f)} \right)}{a_{M}\left( {\theta(f)} \right)}{\psi_{\ell}(f)}{\mathbb{d}f}}}} = {{\sqrt{M}{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{{\psi_{\ell}(f)}}^{2}{\mathbb{d}f}}}} = \sqrt{M}}}} & (36) \end{matrix}$ and interference

$\begin{matrix} {B_{\ell,\ell^{\prime}} = {{\frac{1}{\sqrt{M}}{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{\psi_{\ell}^{*}(f)}{a_{M}^{H}\left( {\theta(f)} \right)}{a_{M}\left( {\theta(f)} \right)}{\psi_{\ell^{\prime}}(f)}{\mathbb{d}f}}}} = {{\sqrt{M}{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{\psi_{\ell}^{*}(f)}{\psi_{\ell^{\prime}}(f)}{\mathbb{d}f}}}} = 0.}}} & (37) \end{matrix}$ Thus, the matched filter introduces no interference regardless of the choice of basis functions, and the SINR is the SNR and is the same for all z_(l)

$\begin{matrix} {{SINR}_{\ell} = {{SNR}_{\ell} = {M\;{\frac{E_{s}}{N_{o}}.}}}} & (38) \end{matrix}$

Plugging equation (35) into equation (30) yields z _(l) =∫g _(l) ^(H)(t) r (t)dt=∫ψ* _(l)(τ)[∫h ^(H)(t−τ) r (t)dt]dτ  (39) so the optimal receiver can be interpreted as a bank of M adjustable delay filters and phase shifters defined by h_(k)(τ) followed by spatial combining and correlation with the basis functions, which corresponds to true time delay beamforming.

A phased array receiver corrects for the phase shift across array 202 in equation (14), but does not perform any temporal equalization to account for the delays, and as a result yields

$\begin{matrix} {{{\underset{\_}{\mathcal{g}}}_{l}(t)} = {\frac{1}{\sqrt{M}}{a_{M}\left( \theta_{O} \right)}{{\Psi_{l}(t)}.}}} & (40) \end{matrix}$ with norm ∥g_(l)∥²=1. This results in signal amplitude

$\begin{matrix} {A_{l} = {{\frac{1}{\sqrt{M}}{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{\Psi_{l}^{*}(f)}{a_{M}\left( \theta_{O} \right)}{a_{M}\left( {\theta(f)} \right)}{\Psi_{l}(f)}\ {\mathbb{d}f}}}} = {\frac{1}{\sqrt{M}}{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{{\Psi_{l}(f)}}^{2}{D_{M}\left( {\theta_{O}\frac{f}{f_{C}}} \right)}\ {\mathbb{d}f}}}}}} & (41) \end{matrix}$ and interference

$\begin{matrix} {B_{l,l^{\prime}} = {\frac{1}{\sqrt{M}}{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{\Psi_{l}^{*}(f)}{a_{M}\left( \theta_{O} \right)}{a_{M}\left( {{{\theta(f)}{\Psi_{l^{\prime}}(f)}{\mathbb{d}f}} = {\frac{1}{\sqrt{M}}{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{\Psi_{l}^{*}(f)}{\Psi_{l^{\prime}}(f)}{D_{M}\left( {\theta_{O}\frac{f}{f_{C}}} \right)}\ {{\mathbb{d}f}.}}}}} \right.}}}}} & (42) \end{matrix}$

In contrast with the matched filter, the SINR of the phased array receiver depends on Δ_(ch) and the choice of basis functions. From FIGS. 5a-5e , it was shown that as Δ_(ch) increases, the value of

$D_{M}\left( {\theta_{o}\frac{f}{f\; c}} \right)$ varies over the bandwidth. From equation (41), this variation results in power loss, and from equation (42), results in interference between the basis functions.

The implementation of the phased array receiver corresponds to removing the filter bank used for space-time equalization in the optimal receiver. In the special case when θ_(o)=i_(o)Δθ, this also corresponds to correlating the i_(o) ^(th) element of the beamspace received signal, r_(b,i) _(o) (t), with the basis functions {ψ_(l)(t)}.

The optimal matched filter receiver indicated by equation (35) in beamspace is

$\begin{matrix} {{{\underset{\_}{\mathcal{g}}}_{b,l}(t)} = {\frac{1}{\sqrt{M}}\left( {{\underset{\_}{h}}_{b}^{*}\Psi_{l}} \right)(t)}} & (43) \end{matrix}$ where the space-time equalization is performed over all M beams. However, as discussed previously, the majority of the channel power is captured by approximately Δ_(ch) beams. Thus, a low-complexity receiver can be designed by processing only these dominant beams. A set

⊂I(M) of dominant beams can be defined as

={iεI(M):σ_(i) ²≧γ}  (44) where σ_(i) ² is a channel power captured by the i^(th) beam from equation (27). A threshold γ is pre-defined and may be defined so that the number of dominant beams are at least equal to the channel spread factor, that is p=|

|≈Δ_(ch), and also capture the majority of the channel power to result in an approximately flat energy density.

Selection of the dominant beams using the threshold and equation (44) is useful in the general case of line-of-sight and/or multipath propagation because it selects all dominant beams up to the point of recovering a desired (and sufficiently large) fraction of the channel power (e.g., 95% of channel power is included so γ is set to capture at least 95% of the total channel power) (setting γ to 95% of the total channel power may result in at most one and possibly no beams being selected. For example, if p=3 as in FIG. 5a-c , each beam captures about 33% of the channel power so no beams may be selected) so that the resulting energy spectral density is nearly flat (and thus, intersymbol interference is minimized) as illustrated in FIGS. 6a and 6b . In essence, each propagation path is mapped to multiple beams (approximately equal to the channel spread factor). Since the channel spread factor is known a priori based on the beam direction, the channel spread factor can also help determine the number of beams associated with each propagation path without using the explicit thresholding of equation (44). In general, due to energy leakage when a fixed beam does not exactly line up with a path direction, the selected number of beams associated with each path is a small multiple, such as 2 to 4, of the channel spread factor.

An optimal receiver's beamspace implementation is shown in FIG. 7. The complex baseband signal r(t) is received at array 202. The complex baseband signal is transformed to the beamspace system model by applying beamforming matrix 306, U_(M), as indicated by equation (7) in the time domain or by equation (8) in the frequency domain. A bank of M filters 700 defined by h _(b,i)(τ) is applied to the signal after the beamspace transformation (U_(M)). The filtered beams may be combined by an additive mixer 702 as understood by a person of skill in the art. The combined beams may be correlated with the chosen basis functions using a correlator 704 to determine the l^(th) test statistic (z_(l)) as understood by a person of skill in the art. The information symbols {s_(l)} transmitted by second transceiver 102 can be detected using the l test statistics (z_(l)), depending on the digital constellation used at the transmitter, as understood by a person of skill in the art. Application of beamforming matrix 306, U_(M), the bank of M filters 700, additive mixer 702, and correlator 704 to the received complex baseband signal r(t) may be performed by processor 302 of receiver 300.

The bank of M filters 700 include a spatial phase shift component and a temporal delay component. For example, the spatial phase shift component may be computed as e^(j2πk(θ) ^(o) ^(−iΔθ)). For example, the temporal delay component may be computed as

$\tau_{k} = {{k\;{\delta\tau}} = {\frac{k\;\theta_{O}}{f_{C}}.}}$ As shown in rig. 7, each filter of the bank of M filters 700 may be defined as

$\frac{1}{\sqrt{M}}{\sum\limits_{k \in {I{(M)}}}{{\mathbb{e}}^{{j\; 2\pi\;{k{({\theta_{O} - {i\;{\Delta\theta}}})}}}\;}W\;\sin\;{{c\left( {W\left( {t + \tau_{k}} \right)} \right)}.}}}$

A near-optimal receiver represented in FIG. 8 can be defined by

$\begin{matrix} {{{\mathcal{g}}_{b,l,i}(t)} = \left\{ \begin{matrix} {\frac{1}{\sqrt{M}}\left( {h_{b;i}^{*}\Psi_{l}} \right)(t)} & {i \in} \\ 0 & {i \notin} \end{matrix} \right.} & (45) \end{matrix}$ with norm

${{{\underset{\_}{\mathcal{g}}}_{b,l}}^{2} = {\frac{1}{M^{2}}{\int{{{\Psi_{l}(f)}}^{2}\left( {\sum\limits_{i \in}{D_{M}^{2}\left( {{\theta(f)} - {i\;{\Delta\theta}}} \right)}} \right){\mathbb{d}f}}}}},$ which results in the following signal amplitude and interference expressions

$\begin{matrix} {A_{l} = {{\frac{1}{\sqrt{M}}{\int{{\Psi_{l}^{*}(f)}\left( {\sum\limits_{i \in}{{H_{b,i}(f)}}^{2}} \right){\Psi_{l}(f)}{\mathbb{d}t}}}} = {\frac{1}{M\frac{3}{2}}{\int_{- \frac{W}{2}}^{\frac{W}{2}}\ {{{\Psi_{l}(f)}}^{2}\left( {\sum\limits_{i \in}{D_{M}^{2}\left( {{\theta(f)} - {i\;{\Delta\theta}}} \right)}} \right){\mathbb{d}f}}}}}} & (46) \\ {B_{l,l^{\prime}} = {\frac{1}{M\frac{3}{2}}{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{\Psi_{l}^{*}(f)}\left( {\sum\limits_{i \in}{D_{M}^{2}\left( {{\theta(f)} - {i\;\Delta\;\theta}} \right)}} \right){\Psi_{l^{\prime}}(f)}\ {{\mathbb{d}f}.}}}}} & (47) \end{matrix}$

When

$p \approx {\Delta_{{ch}^{\prime}}{\sum\limits_{i \in}{\frac{1}{M}{D_{M}^{2}\left( {{\theta(f)} - {i\;\Delta\;\theta}} \right)}}}} \approx M$ over the bandwidth, A_(l)≈√{square root over (M)} and B_(l)≈0. So, in contrast to the phased array receiver, by combining beamforming and space-time equalization, the low-complexity receiver is able to attain near-optimum performance.

As summarized in FIG. 8, the complex baseband signal r(t) is received at array 202. The complex baseband signal is transformed to the beamspace system model by applying beamforming matrix 306, U_(M), as indicated by equation (7) in the time domain or by equation (8) in the frequency domain. The beams p in

are selected using equation (44) and the pre-defined threshold γ. A bank of filters 800 are defined according to equation (45) based on the selected beams p, iε

. The filtered beams (based on the channel responses h_(b,i)(t) as in equation (18)) may be combined by additive mixer 702. The combined beams may be correlated with the chosen basis functions using correlator 704 to determine the l^(th) test statistic (z_(l)) as understood by a person of skill in the art. The information symbols {s_(l)} transmitted by second transceiver 102 can be detected using the l test statistics (z_(l)), depending on the digital constellation used at the transmitter. Application of beamforming matrix 306, U_(M), the bank of filters 800, additive mixer 702, and correlator 704 to the received complex baseband signal r(t) may be performed by processor 302 of receiver 300. The low-complexity near-optimal implementation differs by only performing the filtering and combining over the p beams in

as shown in FIG. 8. Because p≈Δ_(ch), this corresponds to a complexity reduction by a factor of approximately αθ_(o).

Additional functions of a receiver 300 may include downmixing and analog-to-digital (A/D) conversion. As understood by a person of skill in the art, downmixing and A/D conversion may be implemented at different points in the signal processing as illustrated in FIGS. 7 and 8. For example, in an all digital receiver implementation (beamforming matrix 306, the bank of M filters 700 or the bank of filters 800, additive mixer 702, and correlator 704 implemented digitally), downmixing and A/D conversion may be implemented at point “A” indicated in FIGS. 7 and 8, after receipt of the signal by the plurality of antennas 200 of array 202 and before further processing by beamforming matrix 306.

As another example, in a digital receiver with an analog beamforming implementation (beamforming matrix 306 implemented in analog passband, for example, by a lens, and the bank of M filters 700 or the bank of filters 800, additive mixer 702, and correlator 704 implemented digitally), downmixing and A/D conversion may be implemented at point “B” indicated in FIGS. 7 and 8, after processing by beamforming matrix 306 and before further processing by the bank of M filters 700.

As yet another example, in an analog receiver (beamforming matrix 306, the bank of M filters 700 or the bank of filters 800, additive mixer 702, and correlator 704 implemented in analog passband), downmixing and A/D conversion may be implemented at point “C” indicated in FIGS. 7 and 8, after processing by correlator 704 and before any further signal processing.

For illustration, receiver performance can be compared for a sinc pulse single carrier communication system where the N basis functions are delayed versions of the sinc pulse of bandwidth W:

$\begin{matrix} {{{\Psi_{l}(t)} = {\frac{1}{\sqrt{W}}\sin\;{c\left( {W\left( {t - \frac{l}{W}} \right)} \right)}}},{{\Psi_{l}(f)} = {{\mathbb{e}}^{{- j}\; 2\pi\frac{l}{W}f}.}}} & (48) \end{matrix}$ Since |ψ_(l)(f)|²=1, the signal power |A_(l)|² is constant. Similarly, since

${{{\Psi_{l}^{*}(f)}{\Psi_{l^{\prime}}(f)}} = {\mathbb{e}}^{{j2\pi}\frac{l - l^{\prime}}{W}}},$ from equations (42) and (47), the interference terms B_(l,l′) are completely determined by m=l−l′ with B_(m)=B*_(−m). Thus, the interference power is symmetric |B_(m)|²=|B_(−m)|². Results were calculated for an M=61 element,

$\frac{\lambda_{c}}{2}$ spaced array (6 inch array at 60 gigahertz) receiving a signal from φ=55° (θ_(o)=25Δθ). The signal amplitude and interference terms were calculated through numerical integration for the central l=┌N/2┐^(th) pulse. For the largest a considered, for all receivers, the |B_(M)|² values drop below −40 dB of the maximum when |m|≧8. Thus, interference for |m|≧8 can be ignored and the results provide a good assessment of the SINR of any pulse for N>16, except for edge cases (e.g. l=0) where interference will be at most 3 dB lower due to symmetry of the |B_(M)|².

When considering the results, from equation (33) the SINR values at low

$\frac{E_{s}}{N_{o}}$ represent me power loss and from equation (34) the SINR values at high

$\frac{E_{s}}{N_{o}}$ represent the interference caused by pulse distortion. FIG. 9a plots the SINR for the phased array receiver for several values of α(Δ_(ch)), which as expected exhibits severe performance losses as a increases. Even at the relatively low value of α=0.02 (Δ_(ch)=0.5) where there is essentially no power loss, interference causes a 15 dB loss of SINR compared to an optimum value at

$\frac{E_{s}}{N_{o}}$ of 30 dB. Further increasing α only causes more interference, and once α≧0.08 (Δ_(ch)≧2) significant power loss occurs since D_(M) (θ_(o)f/f_(c)) includes nulls in the bandwidth (see FIGS. 5a-e ) resulting in severe power loss.

The low-complexity receiver depicted in FIG. 8 showed much less interference and power loss as shown in FIG. 9b for a 3-beam receiver (p=3) and in FIG. 9c for a 5-beam receiver (p=5). While there is still interference, it was significantly reduced as compared to the phased array receiver. Furthermore, when p≧Δ_(ch), essentially no power loss occurred. FIGS. 9b and 9c also illustrate the ability of the low-complexity receiver depicted in FIG. 8 to deliver near-optimum performance by including a sufficiently large number of beams—at least Δ_(ch). Thus, the low-complexity near-optimal implementation shown in FIG. 8, delivers near-optimum performance with dramatically reduced complexity. The 5-beam receiver outperformed the 3-beam receiver in terms of interference even when Δ_(ch)≦3.

For an OFDM system with N subcarriers the basis functions are

$\begin{matrix} {{{\Psi_{l}(t)} = {\frac{1}{\sqrt{T}}{\mathbb{e}}^{j\; 2\pi\; f_{l}t}}},{0 \leq t \leq {T.}}} & (49) \end{matrix}$ where T=NW and the subcarrier center frequencies are f _(l)=(l−└N/2┘)1/T.  (50)

The Fourier transforms of the ψ_(l)(t) are ψ_(l)(f)=√{square root over (T)}e ^(−πT(f-f) ^(l) ⁾sinc(T(f−f _(l))).  (51)

With the approximation √{square root over (T)}sinc(T(f−f_(l))→δ(f−f_(l)) as T→∞ and from the interference expressions for the phased array and near optimal beamspace receivers, in both cases B_(l,l′)≈0. Thus, for both receivers, the use of OFDM eliminates intersymbol interference. However, the signal amplitudes show a significant difference. From equation (41), the signal amplitude for the phased array receiver is

$A_{l} \approx {\frac{1}{\sqrt{M}}{{D_{M}\left( {\theta_{O}\frac{f_{l}}{f_{c}}} \right)}.}}$ D_(M)(θ_(o)f/f_(c)) occupies approximately 1/Δ_(ch) of the available bandwidth. So only approximately 1/Δ_(ch) of the OFDM subcarriers have a sufficiently high SNR to support communication.

On the other hand, the amplitude for the near optimal B-MIMO receiver is

$A_{l} \approx {\frac{1}{M\frac{3}{2}}{\sum\limits_{i \in {}}{{D_{M}^{2}\left( {{\theta\left( f_{l} \right)} - {i\;{\Delta\theta}}} \right)}.}}}$ When p=|

|≈Δ_(ch),

${\sum\limits_{i \in}{D_{M}^{2}\left( {{\theta\left( f_{l} \right)} - {i\;{\Delta\theta}}} \right)}} \approx {M^{2}.}$ So the near optimal B-MIMO receiver signal amplitude is approximately constant across the subcarriers with A_(l)≈√{square root over (M)}. Thus, in contrast with the phased array receiver, if p≈Δ_(ch) the B-MIMO receiver has a sufficiently high SNR for communication for all the OFDM subcarriers and is able to fully exploit the available bandwidth.

In a system with transmit basis waveforms of the form

$\begin{matrix} {{\Psi_{l}(t)} = {\sum\limits_{n = 0}^{N - 1}{\Psi_{l,n}{u\left( {t - {nT}_{s}} \right)}}}} & (52) \end{matrix}$ where T_(s) is the symbol period,

$W_{s} = \frac{1}{T_{s}}$ is the signaling bandwidth, N≈W_(s)T, ψ_(l)=[ψ_(l,0) . . . ψ_(l,N-1)]^(T), and the {ψ}_(l=0) ^(N-1) are the digital basis waveforms with <ψ_(l),ψ_(l′)>=ψ_(l) ^(H)ψ_(l′)=δ_(l,l′).  (53) The pulse u(t) with Fourier transform U(f)=

{u(t)} represents the pulse shaping filter with bandwidth W≧W_(s) applied to the output of the digital-to-analog converter (DAC) that converts the digitial waveforms to analog waveforms. The inner product of the basis waveforms is

$\begin{matrix} {\left\langle {\Psi_{l},\Psi_{l^{\prime}}} \right\rangle = {{\int_{0}^{T}{{\Psi_{l}^{*}(t)}{\Psi_{l^{\prime}}(t)}\ {\mathbb{d}t}}} = {\sum\limits_{n = 0}^{N - 1}{\sum\limits_{n^{\prime} = 0}^{N - 1}{\Psi_{l,n}^{*}\Psi_{l^{\prime},n^{\prime}}{\int_{0}^{T}{{u^{*}\left( {t - {nT}_{s}} \right)}{u\left( {t - {n^{\prime}T_{s}}} \right)}{\mathbb{d}t}}}}}}}} & (54) \\ {= {\sum\limits_{n = 0}^{N - 1}{\sum\limits_{n^{\prime} = 0}^{N - 1}{\Psi_{l,n}^{*}\Psi_{l^{\prime},n^{\prime}}{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{{U(f)}}{\mathbb{e}}^{{- j}\; 2\pi\frac{n^{\prime} - n}{W_{s}}}\ {\mathbb{d}f}}}}}}} & (55) \end{matrix}$

If the pulse u(t) is chosen to be the root raised cosine filter with bandwidth W=(1+β)W_(s) it can be shown that

$\begin{matrix} {{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{{U(f)}}^{2}{\mathbb{e}}^{{- j}\; 2\pi\frac{n^{\prime} - n}{W_{s}}}\ {\mathbb{d}f}}} = \delta_{n,n^{\prime}}} & (56) \end{matrix}$ and the basis waveforms are orthonormal. Plugging equation (52) into the expression for the matched filter equation (35) yields

$\begin{matrix} {{{\underset{\_}{\mathcal{g}}}_{l}(t)} = {{\frac{1}{\sqrt{M}}\left( {{\underset{\_}{h}}^{*}\Psi_{l}} \right)(t)} = {{\frac{1}{\sqrt{M}}{\int_{0}^{T}{{\Psi_{l}{()}}{\underset{\_}{h}\left( {t -} \right)}\ \mathbb{d}}}} = {\frac{1}{\sqrt{M}}{\sum\limits_{n = 0}^{N - 1}\;{\Psi_{l,n}{\int_{0}^{T}{{u\left( {- {nT}_{s}} \right)}{\underset{\_}{h}\left( {t -} \right)}\ {{\mathbb{d}t}.}}}}}}}}} & (57) \end{matrix}$

From equation (16), the k^(th) element of the integral in (57) can be written as

$\begin{matrix} {{\int_{0}^{T}{{u\left( {- {nT}_{s}} \right)}{{\mathbb{e}}^{{- j}\; 2\pi\; k\;\theta}\  \circ W}\;\sin\;{c\left( {W\left( {t - - {k\;\delta\;}} \right)} \right)}\mathbb{d}}} = {{{\mathbb{e}}^{{- {j2\pi}}\; k\;\theta} \circ {u\left( {t - {k\;\delta} - {n\; T_{s}}} \right)}}.}} & (58) \end{matrix}$

Thus, if the basis waveforms are chosen according to (52) the receiver implementation consists of replacing the sinc filters shown in FIGS. 7 and 8 with u(t), sampling the output of the combiner at t={nT_(s)}_(n=0) ^(N-1), and performing the inner products in the digital domain with the {ψ_(l)}_(l=0) ^(N-1).

Given a multiple-input, single-output (MISO) system with an M_(T)-element ULA transmitter communicating with a single-antenna receiver located at φ_(T,o)

θ_(T,o), the LoS wideband MISO system model (baseband frequency domain) is R(f)= H ^(H)(f) X (f)+W(f); H (f)=a _(M) _(T) (θ_(T)(f))  (59) where R(f) is the received signal, X(f) is the transmitted signal vector, W(f) is noise, and the frequency-dependent beam angle θ_(T)(f) is related to θ_(T,o) as in equation (17). The beamspace LoS wideband MISO system model is

$\begin{matrix} {{{{R(f)} = {{{{\underset{\_}{H}}_{b}^{H}(f)}{{\underset{\_}{X}}_{b}(f)}} + {W(f)}}};{{{\underset{\_}{H}}_{b}(f)} = \left\lbrack {H_{b,i}(f)} \right\rbrack_{i \in {I{(M)}}}}},{{H_{b,i}(f)} = {\frac{1}{\sqrt{M}}{D_{M}\left( {{\theta_{T}(f)} - {i\;{\Delta\theta}_{T}}} \right)}}}} & (60) \end{matrix}$

where X _(b)(f) is the beamspace transmitted signal vector and Δθ_(T)=1/M_(T) is the transmit orthogonal beam spacing.

For MIMO communication between an M_(T) dimensional ULA transmitter and an M_(R) dimensional ULA receiver, the wideband MIMO system model is R (f)=H(f) X (f)+ W (f)  (61) where R(f) is the received signal vector, X(f) is the transmitted signal vector, H(f) is the M_(R)×M_(T) spatial channel frequency response, and W(f) is noise. The beamspace wideband MIMO model is R _(b)(f)=H _(b)(f) X _(b)(f)+ W _(b)(f)  (62) where R _(b)(f) is the beamspace received signal vector, X _(b)(f) is the beamspace transmitted signal vector, H_(b)(f)=U_(M) _(R) ^(H)H(f)U_(M) _(T) is the M_(R)×M_(T) spatial channel frequency response, and W_(b)(f) is beamspace noise. The antenna domain and beamspace LoS wideband MIMO channel models readily follow from the SIMO and MISO models:

$\begin{matrix} {{H(f)} = {{{a_{M_{R}}\left( {\theta_{T}(f)} \right)}{a_{M_{T}}^{H}\left( {\theta_{R}(f)} \right)}{H_{b}(f)}} = {{U_{M_{R}}^{H}{H(f)}U_{M_{T}}} = \left\lbrack {H_{b,i,m}(f)} \right\rbrack_{{i \in {(M_{R})}},{m \in {I{(M_{T})}}}}}}} & (63) \\ {{H_{b,i,m}(f)} = {\frac{1}{\sqrt{M_{R}M_{T}}}{D_{M_{R}}\left( {{\theta_{R}(f)} - {i\;{\Delta\theta}_{R}}} \right)}{D_{M_{T}}\left( {{\theta_{T}(f)} - {m\;{\Delta\theta}_{T}}} \right)}}} & (64) \end{matrix}$ where θ_(T)(f) and θ_(R)(f) are related to θ_(T,o) and θ_(R,o), respectively as in equation (17) and Δθ_(T)=1/M_(T) and Δθ_(R)=1/M_(R) are the orthogonal beam spacings for the transmit and receive ULAs. The LoS path is associated with Δ_(ch,R)=M_(R)αθ_(R) beams at the receiver and Δ_(ch,T)=M_(T)αθ_(T) beams at the transmitter corresponding to a Δ_(ch,R)×Δ_(ch,T) sub-matrix of the beamspace matrix H_(b)(f). The channel matrix for a wideband multipath MIMO channel can be modeled as

$\begin{matrix} {{H(f)} = {\sum\limits_{q = 1}^{N\; p}{\beta_{q}{a_{M_{R}}\left( {\theta_{R,q}(f)} \right)}{a_{M_{T}}^{H}\left( {\theta_{T,q}(f)} \right)}{\mathbb{e}}^{{- {j2\pi}}\; T_{q}f}}}} & (65) \end{matrix}$ where N_(P) denotes the number of paths, and the q-th path is associated with a path gain β_(q), angle of departure φ_(T,q), angle of arrival φ_(R,q), and delay T_(q). The physical angles (φ_(T,q),φ_(R,q))

(θ_(T,q),θ_(R,q)) induce frequency-dependent beam angles (θ_(T,q)(f),θ_(R,q)(f)) as in (17). The beamspace representation of H(f) is given by

$\begin{matrix} {{H_{b}(f)} = {{U_{M_{R}}^{H}{H(f)}U_{M_{T}}} = {{\left\lbrack {H_{b,i,m}(f)} \right\rbrack_{{i \in \;{I{(M_{R})}}},{m \in \;{I{(M_{T})}}}}{H_{b,i,m}(f)}} = {\frac{1}{\sqrt{M_{R}M_{T}}}{\sum\limits_{q = 1}^{N_{p}}{\beta_{q}{D_{M_{R}}\left( {{\theta_{R,q}(f)} - {i\;{\Delta\theta}_{R}}} \right)}{D_{M_{T}}\left( {{\theta_{T,q}(f)} - {m\;{\Delta\theta}_{T}}} \right)}{{\mathbb{e}}^{{- {j2\pi}}\; T_{q}f}.}}}}}}} & (66) \end{matrix}$

Transmit/receive channel dispersion factors can be defined for each path: Δ_(ch,T,q)=M_(T)αθ_(T,q) and Δ_(ch,R,q)=M_(R)αθ_(R,q). Each path is associated with Δ_(ch,R,q) beams at the receiver and Δ_(ch,T,q) beams at the transmitter captured by a corresponding Δ_(ch,R,q)×Δ_(ch,T,q) sub-matrix of the beamspace matrix H_(b)(f).

The MISO antenna domain and beamspace transmit signal vectors X(f) and X _(b) (f) are related to the transmit signal S(f) via X (f)= G _(T)(f)S(f)  (67) X _(b)(f)= G _(T,b)(f)S(f)  (68) where G _(T)(f) and G _(T,b)(f) are the antenna domain and beamspace space-time precoders respectively with ∥G _(T)(f)∥²≦1 and ∥G _(T,b)(f)∥²≦1. For a single information symbol s, the transmit signal is S(f)=sψ(f). The receiver correlates R(f) with ψ(f) to obtain the sufficient statistic for detecting s as

$\begin{matrix} {z = {{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{R(f)}{\Psi^{*}(f)}\ {\mathbb{d}f}}} = {{s{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{E(f)}{{\Psi(f)}}^{2}\ {\mathbb{d}f}}}} + w}}} & (69) \end{matrix}$ where E(f)= H ^(H)(f) G _(T)(f)= H _(b) ^(H),(f) G _(b,T)(f)  (70) is the effective frequency response of the channel and precoder and w˜

(0,No) represents the noise. The optimal, phased array, and near optimal B-MIMO MISO precoders are

$\begin{matrix} {{{\underset{\_}{G}}_{T}^{opt}(f)} = {\frac{1}{\sqrt{M}}{a_{M}\left( {\theta(f)} \right)}}} & (71) \\ {{{\underset{\_}{G}}_{T}^{pa}(f)} = {\frac{1}{\sqrt{M}}{a_{M}\left( \theta_{O} \right)}}} & (72) \\ {{{\underset{\_}{G}}_{T,b}^{no}(f)} = \left\{ \begin{matrix} {\frac{1}{\sqrt{M}}{H_{b;i}(f)}} & {i \in} \\ 0 & {i \notin} \end{matrix} \right.} & (73) \end{matrix}$ where the transmit beam mask

is defined analogously to the receive beam mask

for the near optimal B-MIMO SIMO receiver, which results in effective frequency responses for each precoder

$\begin{matrix} {{E^{opt}(f)} = \sqrt{M}} & (74) \\ {{E^{pa}(f)} = {\frac{1}{\sqrt{M}}{D_{M}\left( {\theta_{O}\frac{f}{f_{c}}} \right)}}} & (75) \\ {{E^{no}(f)}\; = {\frac{1}{M^{3/2}}{\sum\limits_{i \in}{{D_{M}^{2}\left( {{\theta(f)} - {i\;\Delta\;\theta}} \right)}.}}}} & (76) \end{matrix}$

For p=|

|≈Δ_(ch) E^(no)(f)≈M and transmit signals S(f) representing a packet of duration T consisting of N≈TW symbols with energy

modulated onto basis waveforms {ψ_(l)(t)}_(l=0) ^(N-1), the SINR of the sufficient statistics {z_(l)}_(l=0) ^(N-1) is

$\begin{matrix} {{{SIN}\; R_{l}} = {\frac{\frac{E_{s}}{No}{A_{l}}^{2}}{{\frac{E_{s}}{No}{\sum\limits_{l^{\prime} \neq l}{B_{l,l^{\prime}}}^{2}}} + 1}.}} & (77) \end{matrix}$

The signal amplitude Δ_(t) and interference B_(l) terms for a given precoder with an associated effective frequency response are

$\begin{matrix} {A_{l} = {\int_{- \frac{W}{2}}^{\frac{W}{2}}{{E(f)}{{\Psi_{l}(f)}}^{2}{\mathbb{d}f}}}} & (78) \\ {B_{l,l^{\prime}} = {\int_{- \frac{W}{2}}^{\frac{W}{2}}{{E(f)}{\Psi_{l^{\prime}}(f)}{\Psi_{l}^{*}(f)}{{\mathbb{d}f}.}}}} & (79) \end{matrix}$

The processing in LoS MIMO corresponds to simultaneously performing MISO precoding at the transmitter and SIMO combining at the receiver. As with the MISO precoder, the antenna domain and beamspace transmit signal vectors X(f) and X _(b) (f) are related to the transmit signal S(f) using X (f)= G _(T)(f)S(f)  (80) X _(b)(f)= G _(T,b)(f)S(f)  (81) where for a single information symbol S(f)=sψ(f). At the receiver, the output of an antenna domain or beamspace combiner (G _(R)(f) or G _(R,b)(f)) is correlated with ψ(f) to obtain the sufficient statistic for detecting s:

$\begin{matrix} {z = {{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{{\underset{\_}{G}}_{R}^{H}(f)}{\underset{\_}{R}(f)}{\Psi^{*}(f)}{\mathbb{d}f}}} = {{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{{\underset{\_}{G}}_{R,b}^{H}(f)}{{\underset{\_}{R}}_{b}(f)}{\Psi^{*}(f)}{\mathbb{d}f}}} = {{s{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{E(f)}{{\Psi(f)}}^{2}{\mathbb{d}f}}}} + w}}}} & (82) \end{matrix}$ where E(f)= G _(R) ^(H)(f)H(f) G _(T)(f)= G _(R,b) ^(H)(f)H _(b)(f) G _(T,b)(f)  (83) is the effective frequency response of the combiner, precoder, and channel, and w˜

(0,N_(o)∥G_(R)(f)ψ(f)∥²) is the noise. For the optimal, phased array, and near optimal B-MIMO MISO precoders and SIMO combiners, the effective frequency responses are

$\begin{matrix} {\mspace{79mu}{{E^{opt}(f)} = \sqrt{M_{R}M_{T}}}} & (84) \\ {\mspace{79mu}{{E^{pa}(f)} = {\frac{1}{\sqrt{M_{R}M_{T}}}{D_{M_{R}}\left( {\theta_{o,R}\frac{f}{f_{c}}} \right)}{D_{M_{T}}\left( {\theta_{o,T}\frac{f}{f_{c}}} \right)}}}} & (85) \\ {{E^{no}(f)} = {\frac{1}{\left( {M_{R}M_{T}} \right)^{3/2}}{\sum\limits_{i \in}{\sum\limits_{m \in}{{D_{M_{R}}^{2}\left( {{\theta_{R}(f)} - {i\;{\Delta\theta}_{R}}} \right)}{{D_{M_{T}}^{2}\left( {{\theta_{T}(f)} - {m\;{\Delta\theta}_{T}}} \right)}.}}}}}} & (86) \end{matrix}$

The B-MIMO transmit and receive beam masks

and

capture the p_(T)|

|≈Δ_(ch,T) dominant transmit beams and the p_(R)=|

|≈Δ_(ch,T) dominant receive beams where Δ_(ch,T) and Δ_(ch,R) are defined using equation (24). As with the SIMO and MIMO system, for transmit signals S(f) representing a packet of duration T consisting of N≈TW symbols with energy

modulated onto basis waveforms {ψ_(l)(t)}_(l=0) ^(N-1), the SINR of the sufficient statistics {z_(l)}_(l=0) ^(N-1) is

$\begin{matrix} {{{{SIN}\; R_{l}} = {\frac{\frac{E_{s}}{No}{A_{l}}^{2}}{{\frac{E_{s}}{No}{\sum\limits_{l^{\prime} \neq l}{B_{l,l^{\prime}}}^{2}}} + {{{{\underset{\_}{G}}_{R}(f)}{\Psi_{l}(f)}}}^{2}}.}}\;} & (87) \end{matrix}$

The signal amplitude A_(l) and interference B_(l) terms for a given precoder with an associated effective frequency response are

$\begin{matrix} {A_{l} = {\int_{- \frac{W}{2}}^{\frac{W}{2}}{{E(f)}{{\Psi_{l}(f)}}^{2}{\mathbb{d}f}}}} & (88) \\ {B_{l,l^{\prime}} = {\int_{- \frac{W}{2}}^{\frac{W}{2}}{{E(f)}{\Psi_{l^{\prime}}(f)}{\Psi_{l}^{*}(f)}{{\mathbb{d}f}.}}}} & (89) \end{matrix}$

For sparse multipath MIMO systems, the MISO precoding and SIMO combining are applied individually to each of the N_(p) paths. The antenna domain and beamspace transmit signal vectors are

$\begin{matrix} {{\underset{\_}{X}(f)} = {\sum\limits_{v = 1}^{N_{p}}{s_{v}{G_{T,v}(f)}{\Psi(f)}}}} & (90) \\ {{{\underset{\_}{X}}_{b}(f)} = {\sum\limits_{v = 1}^{N_{p}}{s_{v}{G_{T,b,v}(f)}{\Psi(f)}}}} & (91) \end{matrix}$ where G_(T,v)(f) and G_(T,b,v)(f) are the antenna domain and beamspace MISO precoders for the v-th path. The receiver applies N_(p) SIMO antenna domain or beamspace combiners for each path and correlates the output to obtain the sufficient statistics for detecting the {s_(v′)}_(v=1) ^(N) ^(p) as

$\begin{matrix} {z_{v} = {{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{{\underset{\_}{G}}_{R,v}^{H}(f)}{\underset{\_}{R}\ (f)}{\Psi^{*}(f)}{\mathbb{d}f}}} = {{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{{\underset{\_}{G}}_{R,b,v}^{H}(f)}{{\underset{\_}{R}}_{b}(f)}{\Psi^{*}(f)}\ {\mathbb{d}f}}} = {{\sum\limits_{v = 1}^{N_{p}}{s_{v}{\int_{- \frac{W}{2}}^{\frac{W}{2}}{{E_{v,v^{\prime}}(f)}{{\Psi(f)}}^{2}\ {\mathbb{d}f}}}}} + w_{v}}}}} & (92) \end{matrix}$ where E _(v,v′)(f)= G _(R,v) ^(H)(f)H(f) G _(T,v′)(f)= G _(R,b,v) ^(H) H _(b)(f) G _(T,b,v′)(f)  (93) is the effective frequency response of the v-th combiner, v′-th precoder, and channel and w_(v)˜

(0,N_(o)∥G_(R,v′)(f) ψ(f)∥²) represents noise. For the optimal, phased array, and near optimal B-MIMO MISO precoder and SIMO combiner, the effective frequency responses are

E v , v ′ opt ⁡ ( f ) = 1 M T ⁢ M R ⁢ ∑ q = 1 N p ⁢ β q ⁢ D M R ⁡ ( θ R , v ⁡ ( f ) - θ R , q ⁡ ( f ) ) ⁢ D M T ⁡ ( θ T , v ′ ⁡ ( f ) - θ T , q ⁡ ( f ) ) ⁢ ⅇ - j ⁢ ⁢ 2 ⁢ π ⁢ ⁢ q f ( 94 )

E v , v ′ ⁡ ( f ) pa = 1 M T ⁢ M R ⁢ ∑ q = 1 N p ⁢ β q ⁢ D M R ⁡ ( θ R , v - θ R , q ⁡ ( f ) ) ⁢ D M T ⁡ ( θ T , v ′ - θ T , q ⁡ ( f ) ) ⁢ ⅇ - j ⁢ ⁢ 2 ⁢ π ⁢ ⁢ q f ( 95 ) E v , v ′ no ⁡ ( f ) = 1 ( M T ⁢ M R ) 3 2 ⁢ ∑ q = 1 N p ⁢ ∑ i ∈ ⁢ ⁢ ∑ M ∈ ⁢ β q ⁢ D M R ⁡ ( θ R , v ⁡ ( f ) - i ⁢ ⁢ Δ ⁢ ⁢ θ R ) ⁢ D M R ⁡ ( θ R , q ⁡ ( f ) - i ⁢ ⁢ Δθ R ) ( 96 ) D _(M) _(T) (θ_(T,v′)(f)−mΔθ _(T))D _(M) _(T) (θ_(R,q′)(f)−mΔθ _(T))e ^(−j2πτ) ^(q) ^(f).

Here the B-MIMO transmit and receive beam masks are

M_(T) = U_(v^(′) = 1)^(N_(p))M_(T, v^(′)) and

M_(R) = U_(v = 1)^(N_(p))M_(R, v) where

and

are defined for each path as in the LoS channel. In sparse systems, the N_(p) paths are sufficiently separated so that the

{M_(T, v^(′))}_(v^(′) = 1)^(N_(p)) and

{M_(R, v)}_(v = 1)^(N_(p)) are disjoint. So the summation terms in each of the E_(v,v′)(f) (94)-(96) are only approximately nonzero for v=v′=q resulting in E_(v,v′)≈0 and for v′=v and

$\begin{matrix} {\mspace{79mu}{{E_{v,v}^{o}(f)} \approx {\beta_{v}\sqrt{M_{R}M_{T}}{\mathbb{e}}^{{- {j2\pi\tau}_{v}}f}}}} & (97) \\ {\mspace{79mu}{{E_{v,v}^{pa}(f)} \approx {\frac{\beta_{v}}{\sqrt{M_{R}M_{T}}}{D_{M_{R}}\left( {\theta_{o,R}\frac{f}{f_{c}}} \right)}{D_{M_{T}}\left( {\theta_{o,T}\frac{f}{f_{c}}} \right)}{\mathbb{e}}^{{- {j2\pi\tau}_{v}}f}}}} & (98) \\ {{E_{v,v}^{no}(f)} \approx {\frac{\beta_{v}}{\left( {M_{R}M_{T}} \right)^{\frac{3}{2}}}{\sum\limits_{i \in}{\sum\limits_{m \in}{{D_{M_{R}}^{2}\left( {{\theta_{R,v}(f)} - {i\;{\Delta\theta}_{R}}} \right)}{D_{M_{T}}^{2}\left( {{\theta_{T,v}(f)} - {m\;{\Delta\theta}_{T}}} \right)}{{\mathbb{e}}^{{- {j2\pi\tau}_{v}}f}.}}}}}} & (99) \end{matrix}$

Comparing these expressions for E_(v,v)(f) to the wideband LoS MIMO expressions for E(f) (see equation (86)), other than the path gains and the linear phase shift caused by the path delays, wideband multipath MIMO systems with sufficiently separated paths are approximately equivalent to N_(p) LoS MIMO channels with the SINR for each path defined according to equation (87) using the appropriate E_(v,v)(f) to calculate the A_(l) and B_(l,l′).

The nature of space-time processing at the transmitter and/or receiver for single-path LoS (or sparse multipath) SIMO, MISO, and MIMO channels has been defined. As noted earlier, the main feature of wideband high-dimensional systems is that each path gets mapped to multiple beams at the transmitter and/or receiver, and wideband space-time communication occurs over the active set of beams at the transmitter and/or receiver.

For illustration, in a multiuser MIMO system, in which a base station equipped with an M-element array is communicating with K single-antenna mobile user devices over a sufficiently large bandwidth (which may still be less than a 10% fractional bandwidth) so that channel dispersion factor associated with each user is at least one so that wideband processing is used. In practice, communication occurs in two phases: a channel estimation phase, followed by a data communication phase. In the channel estimation phase, the dominant set of beams that need to be processed at the base station are determined. This can be done using sequential transmissions from different mobile user devices, one set at the center (carrier) frequency, one set at the lowest frequency in the band, and one set at the highest frequency in the band. With these measurements, a beamspace mask can be determined for the corresponding mobile user device as the set of beams that carry significant power over the span of frequencies in the bandwidth. A complete beamspace-frequency response matrix can be determined for each mobile user device: that is, the gains and phases for each beam in the mask and over all frequencies (e.g., all tones/carriers in an OFDM system). Once the beamspace-frequency response matrices have been estimated at the base station for each mobile user device, they can be used to determine the beamspace-frequency processing/filtering at the base station for suppressing any interference between the mobile user device signals in the communication phase. For example, a minimum mean squared error criterion can be used for determining the beamspace-frequency precoding at the base station for downlink transmissions from the base station to a specific mobile user device, or for determining the beamspace-frequency processing at the base station for uplink transmissions from the specific mobile user device to the base station.

An illustrative beamforming matrix 306 implemented in analog passband is described in Hybrid Analog-Digital Phased MIMO Transceiver System, U.S. Pat. No. 8,811,511 (filed Sep. 28, 2010), that is assigned to the assignee of the present application and that is incorporated herein by reference in its entirety. Referring to FIG. 10, a transmitter 1000, for example of a transmitting transceiver, includes processor 302, the bank of filters 800, array 202, and beamforming matrix 306, U_(M). The plurality of antennas 200 of array 202 may receive a first signal to compute the columns of beamforming matrix 306, U_(M). Processor 302 may form a signal 1002 that includes one or more information symbols that is sent to the bank of filters 800. The bank of filters 800 may be defined in the frequency domain according to equation (73) based on the selected second plurality of beams p, iε

. For each beam of the selected second plurality of beams as a first beam, a spatial phase shift component and a temporal delay component are computed based on the fixed spatial angle associated with the first beam and the carrier frequency to define the bank of filters 800.

Signal 1002 is filtered through each filter of the bank of filters 800. The filtered signals are provided to the plurality of antennas 200 of array 202 which radiate the filtered signals towards beamforming matrix 306. Beamforming matrix 306 radiates the transmit signal towards the receiver at a receiving transceiver. Beamforming matrix 306 may be implemented using a lens such as a discrete lens array.

The word “illustrative” is used herein to mean serving as an example, instance, or illustration. Any aspect or design described herein as “illustrative” is not necessarily to be construed as preferred or advantageous over other aspects or designs. Further, for the purposes of this disclosure and unless otherwise specified, “a” or “an” means “one or more”. Still further, in the detailed description, using “and” or “or” is intended to include “and/or” unless specifically indicated otherwise. The illustrative embodiments may be implemented as a method, apparatus, or article of manufacture using standard programming and/or engineering techniques to produce software, firmware, hardware, or any combination thereof to control a computer to implement the disclosed embodiments.

The foregoing description of illustrative embodiments of the disclosed subject matter has been presented for purposes of illustration and of description. It is not intended to be exhaustive or to limit the disclosed subject matter to the precise form disclosed, and modifications and variations are possible in light of the above teachings or may be acquired from practice of the disclosed subject matter. The embodiments were chosen and described in order to explain the principles of the disclosed subject matter and as practical applications of the disclosed subject matter to enable one skilled in the art to utilize the disclosed subject matter in various embodiments and with various modifications as suited to the particular use contemplated. 

What is claimed is:
 1. A method of processing a communication signal, the method comprising: receiving a signal by a plurality of antennas; computing, by a processor, a plurality of beam measurement values by transforming the signal received from each antenna of the plurality of antennas to a beamspace representation using a beamforming matrix, wherein the beamforming matrix includes a plurality of columns, wherein each column of the plurality of columns defines a steering vector for a beam that corresponds to a fixed spatial angle relative to a boresight vector extending from the plurality of antennas, wherein the plurality of columns thereby define a plurality of beams; selecting, by the processor, a second plurality of beams from the defined plurality of beams based on a power captured by each beam of the plurality of beams that is derived from the computed plurality of beam measurement values; for each beam of the selected second plurality of beams as a first beam, computing, by the processor, a spatial phase shift component based on the fixed spatial angle defined for the first beam and a carrier frequency; computing, by the processor, a temporal delay component based on the fixed spatial angle defined for the first beam and the carrier frequency; defining, by the processor, a filter using the computed spatial phase shift component and the computed temporal delay component; selecting, by the processor, a beam measurement value associated with the first beam from the computed plurality of beam measurement values; and applying, by the processor, the defined filter to the selected beam measurement value to define a filtered value; and combining the defined filtered values to extract information contained in the received signal.
 2. The method of claim 1, wherein the second plurality of beams includes at least three beams.
 3. The method of claim 1, wherein the second plurality of beams are selected using

={iεI(M):σ_(i) ²≧γ}, where M is a number of the plurality of antennas, σ_(i) ² is a channel power captured by an i^(th) beam of the plurality of beams and γ is a pre-defined threshold.
 4. The method of claim 1, wherein the second plurality of beams are selected based on an angle of a transmitter of the signal relative to the boresight vector.
 5. The method of claim 1, wherein the second plurality of beams are selected as Δ_(ch)=Mα|θ_(o)|, where Δ_(ch) is a number of the second plurality of beams, M is a number of the plurality of antennas, α is a fractional bandwidth of the received signal, and ${\theta_{o} = {\frac{d}{\lambda_{c}}\sin\;\left( \phi_{o} \right)}},$ where d is a spatial distance between adjacent antennas of the plurality of antennas, λ_(c) is a wavelength at the carrier frequency, and φ_(o) is an angle of a transmitter of the signal relative to the boresight vector.
 6. The method of claim 1, wherein the spatial phase shift component is computed using e^(j2πk(θ) ^(o) ^(−iΔθ)), where k is an antenna index value, ${\theta_{o} = {\frac{d}{\lambda_{c}}\sin\;\left( \phi_{o} \right)}},$ where d is a spatial distance between adjacent antennas of the plurality of antennas, λ_(c) is a wavelength at the carrier frequency, φ_(o) is an angle of a transmitter of the signal relative to the boresight vector, i is a beam index value to the first beam, and Δθ=1/M, where M is a number of the plurality of antennas.
 7. The method of claim 6, wherein the temporal delay component is computed as ${\tau_{k} = \frac{k\;\theta_{o}}{f_{c}}},$ where f_(c) is the carrier frequency.
 8. The method of claim 7, wherein the filter is defined using ${\frac{1}{\sqrt{M}}{\sum\limits_{k \in {I{(M)}}}{{\mathbb{e}}^{{j2\pi}\;{k{({\theta_{o} - {i\;{\Delta\theta}}})}}}W\;\sin\;{c\left( {W\left( {t + \tau_{k}} \right)} \right)}}}},$ where W is a bandwidth of the signal and t is a time.
 9. A receiver comprising: a plurality of antennas configured to receive a signal; and a processor configured to compute a plurality of beam measurement values by transforming the signal received from each antenna of the plurality of antennas to a beamspace representation using a beamforming matrix, wherein the beamforming matrix includes a plurality of columns, wherein each column of the plurality of columns defines a steering vector for a beam that corresponds to a fixed spatial angle relative to a boresight vector extending from the plurality of antennas, wherein the plurality of columns thereby define a plurality of beams; select a second plurality of beams from the defined plurality of beams based on a power captured by each beam of the plurality of beams that is derived from the computed plurality of beam measurement values; for each beam of the selected second plurality of beams as a first beam, compute a spatial phase shift component based on the fixed spatial angle defined for the first beam and a carrier frequency; compute a temporal delay component based on the fixed spatial angle defined for the first beam and the carrier frequency; define a filter using the computed spatial phase shift component and the computed temporal delay component; select a beam measurement value associated with the first beam from the computed plurality of beam measurement values; and apply the defined filter to the selected beam measurement value to define a filtered value; and combine the defined filtered values to extract information contained in the received signal.
 10. The receiver of claim 9, wherein the second plurality of beams includes at least three beams.
 11. The receiver of claim 9, wherein the second plurality of beams are selected using

={iεI(M):σ_(i) ²≧γ}, where M is a number of the plurality of antennas, σ_(i) ² is a channel power captured by an i^(th) beam of the plurality of beams and γ is a pre-defined threshold.
 12. The receiver of claim 9, wherein the second plurality of beams are selected based on an angle of a transmitter of the signal relative to the boresight vector.
 13. The receiver of claim 9, wherein the second plurality of beams are selected as Δ_(ch)=Mα|θ_(o)|, where Δ_(ch) is a number of the second plurality of beams, M is a number of the plurality of antennas, α is a fractional bandwidth of the received signal, and ${\theta_{o} = {\frac{d}{\lambda_{c}}{\sin\left( \phi_{o} \right)}}},$ where d is a spatial distance between adjacent antennas of the plurality of antennas, λ_(c) is a wavelength at the carrier frequency, and φ_(o) is an angle of a transmitter of the signal relative to the boresight vector.
 14. The receiver of claim 9, wherein the spatial phase shift component is computed using e^(j2πk(θ) ^(o) ^(−iΔθ)), where k is an antenna index value, ${\theta_{o} = {\frac{d}{\lambda_{c}}{\sin\left( \phi_{o} \right)}}},$ where d is a spatial distance between adjacent antennas of the plurality of antennas, λ_(c) is a wavelength at the carrier frequency, φ_(o) is an angle of a transmitter of the signal relative to the boresight vector, i is a beam index value to the first beam, and Δθ=1/M, where M is a number of the plurality of antennas.
 15. The receiver of claim 14, wherein the temporal delay component is computed as ${\tau_{k} = \frac{k\;\theta_{o}}{f_{c}}},$ where f_(c) is the carrier frequency.
 16. The receiver of claim 15, wherein the filter is defined using ${\frac{1}{\sqrt{M}}{\sum\limits_{k \in {I{(M)}}}{{\mathbb{e}}^{{j2\pi}\;{k{({\theta_{o} - {i\;{\Delta\theta}}})}}}W\;\sin\;{c\left( {W\left( {t + \tau_{k}} \right)} \right)}}}},$ where W is a bandwidth of the signal and t is a time.
 17. A transmitter comprising: a plurality of antennas configured to receive a first signal and to transmit a second signal; and a processor configured to compute a plurality of beam measurement values by transforming the first signal to a beamspace representation using a beamforming matrix, wherein the beamforming matrix includes a plurality of columns, wherein each column of the plurality of columns defines a steering vector for a beam that corresponds to a fixed spatial angle relative to a boresight vector extending from the plurality of antennas, wherein the plurality of columns thereby define a plurality of beams; select a second plurality of beams from the defined plurality of beams based on a power captured by each beam of the plurality of beams that is derived from the computed plurality of beam measurement values; for each beam of the selected second plurality of beams as a first beam, compute a spatial phase shift component based on the fixed spatial angle defined for the first beam and a carrier frequency; compute a temporal delay component based on the fixed spatial angle defined for the first beam and the carrier frequency; and define a filter using the computed spatial phase shift component and the computed temporal delay component; filter a transmit signal through each of the defined filters; and provide the filtered transmit signal to the beamforming matrix as part of forming the second signal.
 18. The transmitter of claim 17, wherein the second plurality of beams are selected as Δ_(ch)=Mα|θ_(o)|, where Δ_(ch) is a number of the second plurality of beams, M is a number of the plurality of antennas, α is a fractional bandwidth of the received signal, and ${\theta_{o} = {\frac{d}{\lambda_{c}}{\sin\left( \phi_{o} \right)}}},$ where d is a spatial distance between adjacent antennas of the plurality of antennas, λ_(c) is a wavelength at the carrier frequency, and φ_(o) is an angle of a transmitter of the signal relative to the boresight vector.
 19. The transmitter of claim 17, wherein the second plurality of beams are selected using

={iεI(M):σ_(i) ²≧γ}, where M is a number of the plurality of antennas, σ_(i) ² is a channel power captured by an i^(th) beam of the plurality of beams and γ is a pre-defined threshold.
 20. The transmitter of claim 17, wherein the beamforming matrix is implemented as a lens, wherein the filtered signal is provided to the lens by radiating the filtered signal from the plurality of antennas. 